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Abstract 

The mixing of n — rf or n — n' — G is of a great theoretical interest, because it concerns many 
aspects of the underlying dynamics and hadronic structure of pseudoscalar mesons and glueball. 
Determining the mixing parameters by fitting data is by no means trivial. In order to extract the 
mixing parameters from the available processes where hadrons are involved, theoretical evaluation 
of hadronic matrix elements is necessary. Therefore model-dependence is somehow unavoidable. 
In fact, it is impossible to extract the mixing angle from a unique experiment because the model 
parameters must be obtained by fitting other experiments. Recently BR(D — > n + I + fi) and 
BR{D S — > 77(77') +Z + i>i) have been measured, thus we are able to determine the n — ??' mixing solely 
from the semileptonic decays of D-mesons where contamination from the final state interactions 
is absent. Thus we hope that the model-dependence of the extraction can be somehow alleviated. 
Once BR{D — > rf + I + v\) is measured, we can further determine all the mixing parameters for 
rj — n' — G. As more data are accumulated, the determination will be more accurate. In this work, 
we obtain the transition matrix elements of D( s ) — > rf'' using the light-front quark model whose 
feasibility and reasonability for such processes have been tested. 
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I. INTRODUCTION 



The mixing among pseudoscalar mesons and glueballs is of great theoretical interests and 
significance for understanding the dynamics and hadronic structures. Study on the mixing 
has lasted for several decades, not only because of its importance, but also the difficulties 
caused by both theoretical and experimental aspects. As is well understood, the mixing is 
caused by the QCD anomaly and related to the chiral symmetry breaking [l, Q|. Definitely, 
one would be able to gain a better insight into the dynamics, if the mixing parameters 
are more accurately determined. Many measurements on the processes where rj and rj' are 
involved, have been carried out to fix the mixing parameters. The mixing of 77— 7/ is described 
in different forms, i.e. with different eigen-bases. 

In the SU(3) quark model 77 and rj are the mixtures of r] and 77s which are SU(3) singlet 
and octet states respectively [3j-l5(, 
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and the mixing angle 9 was fitted in a range of —10° to — 23 [4j. Bhattacharya and Rosner's 
recent work [6| suggests that the data of the decays of D°, D ± and D s into two light 
pseudoscalar mesons seem to favor a smaller octet-singlet mixing angle —11.7°. 

In the quark content basis, r/o and rj$ can be further written as mixtures of rj q and r] s , 






(2) 

where i] q = -^(uu + dd) and r] s = ss. In Refs. |7j-l9|, the mixing of 77 and 7/ is written as 

(3) 

where cos <b = \ \ cos 9 — J% sin 9 and sin 6 — \ § cos 9 + \ | sin 9. We refer this mixing as 
scenario-I in this work. The mixing angle 9 = —11.7° [6| corresponds to (f> = 43.0° in this 
scenario. 

In the history, there have been various ways to determine the mixing angle (s). Feldmann 
et al. [k]] summarized the issue and listed several possibilities to extract the mixing angle(s) 
from the experimental data. By fitting the ratio of nonleptonic decay widths r^[j/v>^pf 




the authors of Ref. [U( obtained = (39.9 ± 2.9)°, while by fitting the decay widths of 



T(r]' — > P7) and T(p — > 777), the value of was set as cf> = (35.3 ± 5.5)°. In terms of 
r[j/v^rpj] ^ = (43-1 ± 3.0)° was obtained. Using the ratio of decay widths of p^"^ e e ^ fl2| 
and the form factors given in Ref. 13| <j) is fixed as <fi = (41.3 ± 5.3)° 10] . With the cross 
section of scattering processes n~p —> r/'n and n~p — > r/n one had (J) = (36.5 ± 1.4)° [yj 
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and = (39.3 ± 1.2)° [15]. The Crystal Barrel Collaboration 16] measured the ratio of 
the annihilation processes pp —> 7/+meson(7r°, rj or u) and pp — > ?7+meson(7r , rj or u) and 
achieved 4> = (37.4 ± 1.8)°. Based on the ratio ^fjj^z^^ the value of (p was obtained as 
(39.0 ± 1.6)° 10j. The results seem a bit dispersive, but they are consistent with each other 
for the accuracy the present experiments can reach. Feldmann et al. obtained a weighted 
average of the 4> value as (39.3 ± 1.0)°. Moreover, is not the end of the story, since the QCD 
anomaly, if it causes the mixing between 77 and 7/, also results in a mixing of 77, 7/ with 
gluonium of the quantum number h . 



As one extends the picture to invo 



scenario-II, was suggested in Refs. 
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where |q) = | gluonium) is a pure gluonium state and the physical state G was identified as 



77(1405) [l9j. It is noted that the definition of the mixing angles in Ref. 19| is different from 
that adopted in Ref. j2o|. In Ref. {22} the x 2 scheme was used to obtain 0g ; 0' an d other 
parameters by fitting the data of several radiative decays such as u — » 777, p — > 777 23[ and 
oj — > 7r°7. The results are sin 2 0G = 0.115 ± 0.036 and <p' = (40.4 ± 0.6)° if one sets sin 2 0c as 
a free parameter. 

Recently, the branching ratios of D s — > rj{rj')e + v e and D + — > r]e + v e have been measured 
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r]'e + v e ). Even 



25] and the collaboration 25| obtains a new upper limit of B1Z(D 
though the missing energy of the produced neutrino may bring up certain uncertainties, 
semi-leptonic decay modes have an obvious advantage over the radiative and non-leptonic 
decays. Since the contamination from the final state interaction is absent in semi-leptonic 
decays, the theoretical calculation is more reliable and the physical quantities extracted from 
data, such as the mixing angles, would be closer to reality. This is an ideal opportunity to 
determine the structure of 77(77'), e.g. extract the mixing angle(s) of 77 — 77' (77 — 77' — G) from 
data. This advantage motivates us to study the mixing via the semi-leptonic decays alone. 

As was indicated above, in most cases, the mixing angle and the model parameters were 
not simultaneously determined by fitting solely one type of data. Since all the relevant 
processes involve hadron transitions, one needs to evaluate the hadronic transition matrix 
elements which are fully governed by the non-perturbative QCD, in terms of phenomeno- 
logical models. Therefore extraction of the mixing is somehow model dependent. However, 
if one can determine the mixing angle (s) solely from one type of processes, the model de- 
pendence would be alleviated because he can just employ one phenomenological model to 
deal with all the concerned reactions, e.g. the model parameters and the mixing angles are 
fixed altogether. Indeed, in this way, we can expect that the model dependence which is 
unavoidable, could be reduced to minimum. 
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Our strategy is following. We will concentrate on the study of the semi-leptonic decays of 
D and D s . At the tree- level, only the dd( ss) component of 77, 77' contributes to the transition 
D + — > 77(77') ( D s — > 77(77')). The amplitude at the quark level can be obtained in terms of the 
weak effective theory, so the key point is to calculate the hadronic transition matrix elements. 
Concretely, we are going to use the light front quark model (LFQM) to evaluate the hadronic 
transition matrix elements. It is believed that the LFQM is a relativistic model which has 
obvious advantages for dealing with hadronic transitions where light hadrons are involved 



26, 



271 ] . The light-front wave function is manifestly Lorentz invariant and expressed in terms 



of the fractions of internal momenta of the constituents which are independent of the total 
hadron momentum. Applications of this approach have been discussed in some details by 



the authors of Ref. [29| and their results are in good agreement with data. Moreover, in 
Ref. [3^] we explored the structure of /o(980) via the transition D s — > / (980) in the LFQM. 
However, in the approach, there are a few parameters to be fixed: the /3 values in the hadron 
wave-functions. 

In our scheme, we let the mixing angle and /3 be free parameters, when we fit the data, 
we obtain the mixing angle and the /3 values simultaneously. In this way, we "extract" the 
mixing angle directly from the data on the semileptonic decays of D and D s . At the present 
BTZ(D + — > rje + u e ) and BTZ(Df — > r)(r]')e + v e ) have been measured, but the data are not 
enough to determine mixing angles of 77, 77' and glueball. Once B1Z{D + — > r]'e + b , e ) is measured 
in the future at BES, it would be optimistic that one will be able to further determine the 
mixing structure which is extremely valuable for getting insight into the physics picture of 
light hadrons and probably the glueballs. Thus so far, we only concern the scheme I where 
only 77 — 77' mixing exists. 

This paper is organized as follows: after the introduction, in section II we present the form 
factors of D + — > 77(77') (D s — >■ 77(77')) which are evaluated in the LFQM, then we calculate the 
corresponding decay rates and by fitting the data we gain the mixing angle and the model 
parameters altogether. Section III is devoted to our conclusion and discussions. 



II. CALCULATION OF THE WIDTHS OF D + r](rj')e + u e AND D s -> r)(rf')e+u e IN 
LFQM 

In this section we are going to calculate the decay widths of D + — > ri(r]')e + i> e and D s — > 
rj{rj')e + v e in terms of the LFQM. The crucial task is to evaluate the form factors of D + — > 



r](r]')(D s — y 77(77')). In Ref. 31[, we studied D s — > 77(77') in the LFQM, thus the corresponding 



formulas can be directly used in this work. The transition diagram is shown in Fig. [TJ 
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FIG. 1: Feynman diagram for meson transition amplitude 



A. Formulations 



For being self-content, we list some key formulaes given in Refs. 26|, |27| here. The form 
factors for P — > P transition are defined as 



(P(P")\A,\P(P')) = f + (q 2 )V, + Uq 2 )q». 

It is convenient to redefine them as 

/ M 12 - M" 2 \ M' 2 - M" 2 

(P(P»)\AJP(P>)) = V, - M - ^ < 



where q = P' — P" and V = P' — P" . The relations among the quantities are 



Fi{q 2 ) = U(q 2 ), F (q 2 ) = f + {q 2 



q-V 



(5) 



(6) 



(7) 



Functions f±iq 2 ) can be calculated in the LFQM and their explicit expressions were 
presented as |27| . 



U{q 2 ) 



Uq 2 



Xl {M' 2 + M{ 2 ) - x 2 q 2 + x 2 (m[ - mf ) 

2" 



16tt 3 J - ^x 2 N' l N'{ 
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N c r . .„ . 2hLh" 
16tt 3 J ~ ^x 2 N{N[ 

- 2 q — (p' 2 + 2 ip/± ' q±)2 ) - 2 {p>± ' q±)2 + W± ' q±) 



f dx 2 d 2 p' ± F x {x\X 2 M 12 + p' ± + m' l m 2 + {—m'[ + m 2 )(x 2 m' 1 + 

M" 2 - x 2 {q 2 + q ■ V) 



x 1 m 2 ) 



q \ q ) Q Q 

-(x 2 - x 1 )M' 2 +2 Xl M 12 - 2(m[ - m 2 )(m / 1 + m'()] } . (8) 

where m' 1; m" and m 2 are the corresponding quark masses, M' and M" are the masses of 
the initial and final mesons respectively. All other notations can be found in the Appendix. 
The form factors are parameterized in a three-parameter form as 

F(0) 



F(q 2 



(9) 



5 



where F(q 2 ) represents the form factors F±, Fq, and F(0) is the form factor at q 2 = 0; M is 
the mass of the initial meson. The three parameters F(0) and a, b are fixed by performing 
a three-parameter fit to the form factors which are calculated in the space-like region and 
then extended to the physical time-like region. 

For semileptonic decay of a pseudoscalar meson (D or D s ) into a pseudoscalar meson, i.e. 
P(P') -> P(P")lv, the differential width is 



where q = P' — P" is the momentum transfer and q 2 is the invariant mass of the lepton- 
neutrino pair; p is the final meson momentum in the D or D s rest frame and 



(M 2 - (M f - V¥) 2 ) (m 2 - {M f + y/q*) 

p=\p»\ = ±i . ^ . — . (lr 

Mf denotes the mass of the produced meson. It is noted that the differential width is 
governed by only one form factor Fi(q 2 ) because we neglect the light lepton masses. 



B. Calculating the decay rates of D + (D S ) — > r)(r/)e + i> e and determining the mixing 
angle 

Now, let us extract the mixing angle in scenario-I by means of the experimental data of 
BU(D + r]e + u e ) = (13.3±2.0±0.6) x KT 4 , BK(D a -> r]e + u e ) = (2.48±0.29±0.13) x 10~ 2 
and BK(D S rfe + v e ) = (0.91 ± 0.33 ± 0.05) x 10~ 2 . 

First we need to calculate the form factors F® ^ Ds ^(q 2 ) and F® ( ' Ds ^' 1 (q 2 ) using the 
formulas given in section III Al Since d quark in D + appears in the final meson, only the 
dd component of r)(j]') contributes to the transition of D + — > r](ri')lui i.e. the mixing angle 
is included in the form factors. The similar situation for the transition D s — > r)(rf)lvi was 
discussed in 3J| but we let the mixing angle be free. 

In this work, the input parameters for quark masses are directly taken from 27] as m u = 



0.26 GeV, m s = 0.37 GeV, m c = 1.4 GeV. From the experimental results for decay constants 
f% p = 0.221 MeV, fp p = 0.27 MeV, the parameters for (3 Da and (3 D are fixed to be 



(5 Ds = 0.592 GeV, (5 D = 0.499 GeV [31]. 



Totally, there are five free parameters (3%, (3^,, (3L (3Z and to be fixed. It seems that we 
do not have enough equations to determinate all these parameters. However, the (3p OV s is a 
parameter in the Gaussian wave function whose value only depends on the quark flavor (q 
or s), but not on the characters of the hadron, therefore two relations (3^ = (3^, and (3^ = (3^, 
hold. The relations are also deduced in the Appendix B where we do not make any special 
assumption on them. By these relations we reduce the number of unknowns to three. From 
the three decay modes, we obtain the values of (f> and (3^,-j, The mixing angle is fit 
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FIG. 2: The dependence of the form factor F\(q 2 ) on q 2 at <fi = 39.9°. (a) D + — > r/ (the solid line) 
and D s — > -q (the dotted line); (b) D + — > rj (the solid line) and D s — > rj' (the dotted line) 

to be = (39.9 ± 2.6)° where the errors are from the experiment. The parameters /3 are 
% {rf) = 0.398 GeV, P' {rf) = 0.453 GeV. 

Now, let us discuss the theoretical uncertainties in our used light front quark model. Here, 
we only address the errors coming from the phenomenological parameters. The numerical 



results are not sensitive to the variance of /^S^n and fi^y In [SjJ], it is found that the 
variations of decay constants of fo and fo s with parameters is nearly a linear relation. 
One can ascribe the uncertainties from /?£>, /3d s to experimental errors of fr>, fr> s . For the 
quark masses, we made a variation of their values by 10% and explore the uncertainties of 
predictions. We found that the numerical numbers are not sensitive to m^, m s . In particular, 
with variations of m^, m s , the errors of (f) is only 0.7°. The uncertainties from charm quark 
mass m c is in a controllable region, A0 = 2.1°. Combing them, we have an error 2.3° in 
determination of 0, this is the main error of our approach. After considering the experimental 
error, the mixing angle is determined to be = (39.9 ± 2.6(exp) ± 2.3(the))°. 

SU(3) breaking is a substantial effect in rj — rj' mixing. In our approach, we have included 
their effects in differences of quark masses and m s ; (3 parameter differences in /^j and j3^. 
Using these values, we plot the dependence of the form factor F^q 2 ) on q 2 in the physical 
region of q 2 > in Fig. 2 and we estimate B1Z(D + — > r]'e + u e ) = (2.12 ± 0.23(exp) ± 
0.20(the)) x 10~ 4 which is lower than the experimental upper bound BTZ(D + — > rj'e + v e ) < 



3.5 x 10' 



25]. 



Using the formula of Eq. ( 1B1|) in Appendix and the above parameters we calculate the 
central value of the decay constants f-> = 80 MeV, = 113 MeV, fi = 67 MeV and 
ft = 145 MeV. These values are obtained by fit to the experiment. They are very close 
to the results /« = 78 MeV, f° = 113 MeV, = 64 MeV and f°, = 141 MeV given in 
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Ref. 32J where two mixing angles 6$ = —9.1° and 6$ = —22.2° are used for the mixing of 77 
and 77'. It is noted that our result extracted solely from the data on the semi-leptonic decays 
of D and D s is consistent with the previously determined weighted value 39.3° 10| and 41.4° 



22| | within a reasonable error tolerance. In Refs. lOj, [13] the authors also used the data 
T(D S — > r)'ev)/T(D s — > r\ev) to determine the mixing angle. In this work, we employ an 
alternative model to evaluate the hadronic transition matrix elements. Moreover, nowadays 
there are more data available, which enable us to make more accurate estimation on the 
77 — 77' mixing. 

In principle we will be able to fix the mixing angles 0' and 0g in Eq.(jl]) simultaneously 
as long as a sufficiently large database on the semileptonic decay modes is available. 



III. CONCLUSION 

Because of absence of the final state interactions, the semileptonic decays have obvious 
advantages for determining the properties of the produced light hadrons, such as the structure 
of /o(980), 77 — 77' mixing and even a mixing of pseudoscalar mesons with glueball, over other 
modes from the theoretical aspect. Indeed, moreover, for the semileptonic decays, one only 
needs to calculate the hadronic transition elements where only one final hadron is involved. 
Instead, besides the complication caused by the final state interaction, the evaluation of 
hadronic matrix elements of the non-leptonic decays is much more difficult because there are 
two hadrons in the final state. To do the job, the factorization scheme may be invoked and 
more ambiguities are raised by the scheme. 

In this work we extract the mixing angle of 77 and 77' by the semileptonic decays D — > r/e + z/ e 
and D s — > ri(r]')e + i> e . At first we calculate the form factors of D — > rj and D s — > 77(77') in 
the LFQM where the mixing angle of 77 and 77' and the parameters /3s were free. Then we 
compute the rates of the semileptonic decays D — > r]e + u e and D s — >■ r](r]')e + iy e . Equating our 
theoretical derivations with the data we obtain = (39.9 ± 2.6(exp) ± 2.3(the))°, ^jLn = 
0.398 GeV and 0° {r)l) = 0.453 GeV. We estimate BTZ(D+ r]'e+u e ) = (2.12 ± 0.23(exp) ± 
0.20(the)) x 10~ 4 . 

As indicated above, there are three free parameters including the 77 — 77' mixing in our 
model-dependent calculations, and we fix them simultaneously by the data of BR(D — > 
r]e + u e ) and BR(D S — > ri(r]')e + i> e ) in our scheme I. If in the future B1Z(D + — > r]'e + u e ) can 
be accurately measured (not just an upper bound), we will be able to test our theoretical 
prediction of (2.12±0.23(exp)±0.20(the)) x 10~ 4 which is made in our scheme I. If the data are 
consistent with our estimation, it implies that our scheme I is valid, e.g the mixing between 
77—77' with gluonium states is small. Instead, if there is an observable discrepancy between our 
theoretical prediction on B1Z(D + — > r]'e + i> e ) and the data, one needs to further consider the 
scheme II where the mixing between 77 — 77' is accounted, and by the data, we can determine 
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the mixing angles and 0g in eq.(4). We wish our experimental colleagues to carry out 
measurements on D + — > 7]'e + v e and further improve the accuracy of the measurements of 
semileptonic decays. 

Of course, the branching ratios of semileptonic decay modes are smaller than that of non- 
leptonic modes, and due to the existence of a neutrino (antineutrino), the missing energy 
would make the event re-construction more difficult. However, the obvious advantage of the 
semileptonic decays for determining the properties of the light hadrons, would be the reason 
to urge our experimental colleagues to conduct more precise measurements. It is worthwhile. 
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Appendix A: Notations 

Here we list some variables appearing in the context. The incoming (outgoing) meson in 
Fig. [1] has the momentum P'^"' = p'^"^ + P2 where p'^"^ and P2 are the momenta of the 
off-shell quark and antiquark and 



p'+ = Xl P ,+ , p+ = x 2 P' + , 

Pl± = X l P ± + P±> P2± = %2P'± ~ P±, 

with Xi and p' ± are internal variables and X\ + — 1. 
The variables Mq, Mq, h' p and N[ are defined as 



(Al) 



Mi 



12 




+ 



p'l + ml 




o — 



(A2) 




(A3) 



where 




(A4) 



with p' z 



2 



2x 2 ML ' 



\ =x l (M' 2 -M' Q 2 ). 



(A5) 
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Appendix B: the relations between and /S^n 

For a pseudoscalar meson the decay constant can be evaluated 

f p = / dx 2d 2 p'± ~ 4(7^2 + m 2 xi), (Bl) 



so the parameters /^ ( , } , /S^V, can be determined by the decay constants of /^/) and /^/) 
which are defined in Ref. |10|] 

= /, cos (j) , /" = sin , 
/J = /,sin0, # = / s cos0. (B2) 

where g represents u or d quark. 

Since 77 = r] q cos 4> — rj s sin and 7/ = r/ g sin <p + r] s cos 0, we get the relations (3^ = 
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